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Holonomic gates for quantum computation are commonly considered to be robust against certain kinds of 
parametric noise, the very motivation of this robustness being the geometric character of the transformation 
achieved in the adiabatic limit. On the other hand, the effects of decoherence are expected to become more 
and more relevant when the adiabatic limit is approached. Starting from the system described by Florio et al. 
[Phys. Rev. A 73, 022327 (2006)], here we discuss the behavior of non ideal holonomic gates at finite operational 
time, i.e., far before the adiabatic limit is reached. We have considered several models of parametric noise and 
studied the robustness of finite time gates. The obtained results suggest that the finite time gates present some 
effects of cancellation of the perturbations introduced by the noise which mimic the geometrical cancellation 
effect of standard holonomic gates. Nevertheless, a careful analysis of the results leads to the conclusion that 
these effects are related to a dynamical instead of geometrical feature. 

PACS numbers: 03.67.Lx, 03.65.Vf, 03.65.Yz 
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I. INTRODUCTION 

One of the most important challenges for the realization 
of quantum information tasks is the implementation of quan- 
tum logicgates that are robust against unwanted perturba- 
tions Jll 01 ■ Two kinds of perturbations with qualitatively 
different features can be distinguished: the first kind has a 
purely quantum nature, and it is induced by the interaction 
of the quantum system implementing the logic gate with the 
environment; the second kind has instead a classical nature, 
and it is caused by the presence of instrumental noise in the 
'external parameters' used to control the system. The un- 
wanted interaction with the environment is the source of the 
phenomenon known as quantum decoherence |3|]. The effects 
of this interaction can be modeled by means of suitable 'mas- 
ter equations' (i.e. evolution equations) for the density matrix 
of the quantum system implementing the logic gate; at least 
in the Markovian regime 113511 . they are negligibly small if the 
operational time of the logic gate is short enough. The classi- 
cal perturbations stem from an unavoidable noisy component 
intrinsic in the external driving fields (e.g. laser beams Jit]) 
that can be usually regarded as classical fields; hence, it is es- 
sentially due to instrumental instability. The effects of these 
perturbations can be evaluated by studying standard (non- 
autonomous) Schrodinger equations where the instrumental 
noise is taken into account by suitably modeling the noisy 
components of the classical parameters (e.g. the field ampli- 
tude) associated with the external driving fields. 

Among the several strategies for realizing quantum logic 
gates discussed in the literature, a prominent position is held 
by holonomic gates. They were first proposed by Zanardi 
and Rasetti Jfjl (see also Ref. (0]), and rely on the theory of 
holonomy and of the associated holonomy groups in princi- 
pal fiber bundles JH], a subject which is familiar to theoretical 
physicists due to the central role played in gauge theories ^ 
and in the well-known phenomenon of abelian ifioll and non- 
abelian 1 1 1'] adiabatic phases. Actually, a holonomic gate can 
be regarded as a straightforward application of the theory of 



non-abelian adiabatic phases to quantum computation. 

Since the very beginning, holonomic gates were considered 
to be intrinsically robust against classical noise 11211 . thanks to 
the geometric features of holonomy in Hilbert bundles. As we 
will briefly recall below, three main ingredients are needed in 
order to realize such holonomic gates. 

The first ingredient is a suitable physical system described by 
a quantum Hamiltonian depending on some set of parameters, 
these parameters being associated with the external (classical) 
driving fields that are assumed to be experimentally control- 
lable functions of time; the unavoidable instrumental instabil- 
ity (stochastic noise) affecting the driving fields is the source 
of the classical noise — we will call it parametric noise, in 
the following — that has been mentioned above. 
The second ingredient consists in selecting a suitable 
eigenspace of the given Hamiltonian — an eigenspace de- 
pending smoothly on the external parameters, hence actu- 
ally an iso-degenerate family of eigenspaces; let us call them 
the family of relevant eigenspaces — and in fixing in the 
parameter space an 'initial point' and a loop through this 
point. To such a loop corresponds an excursion of the 
parameter-dependent Hamiltonian (hence, of its eigenprojec- 
tors) and a certain ideal unitary transformation in the encod- 
ing eigenspace, namely, that particular relevant eigenspace 
fixed by the initial (and final) point of the loop in the param- 
eter space. This ideal transformation is determined by Kato's 
adiabatic evolutor associated with the given Hamiltonian and 
with the chosen loop in the parameter space, and it has a sim- 
ple geometric interpretation as a holonomy phenomenon (geo- 
metric phase). The ideal unitary transformation plays a central 
role in Kato's formulation of the adiabatic theorem lfl3ll ap- 
plied to our context. Indeed, the external parameters are con- 
trollable functions of time and in the adiabatic limit — i.e., in 
the limit where the loop in the parameter space is covered in a 
operational time tending to infinity — the real evolution over 
the operational time, determined by the given physical Hamil- 
tonian, becomes cyclic in the encoding eigenspace and, apart 
from an irrelevant overall 'dynamical phase factor', coalesces 
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in this subspace with the ideal unitary transformation. We 
stress that the ideal unitary transformation should be thought, 
in our context, as an ideal quantum gate whose behavior can 
be, in general, only approached by a non-ideal quantum gate 
corresponding to the real evolution over a suitably large, but 
finite, operational time. 

Accordingly, the third ingredient is the choice of a suitable op- 
erational time — which will be called balanced working time, 
in the following — for the real quantum gate. This time span 
must be short enough to achieve a fast quantum computer and 
to avoid the ravages of decoherence, but long enough to justify 
the adiabatic approximation (i.e. to approach the behavior of 
the ideal quantum gate) which is at the root of the appearing 
of geometric phases Dal . Hence, a balanced working time is 
determined by a touchy trade-off between two competing and 
not necessarily compatible demands. 

The problem of robustness of holonomic gates against para- 
metric noise has been studied in both the abelian and non- 
abelian case with different approaches Ennui. In these 
papers, the effects of random perturbations of the control pa- 
rameters are considered. It is worth noticing, however, that 
such effects are evaluated with the adiabatic limit already be- 
ing performed, thus essentially confirming quantitatively the 
standard qualitative geometric argument usually adopted to 
support the robustness of holonomic gates, argument which 
will be recalled later on. 

As holonomic gates are generally considered to be a priori 
robust against parametric noise, attention has mainly focused 
on the study of decoherence effects lfl9ll20ll2lll22Tl and on the 
possibility of partially suppressing them 112311 . These investi- 
gations show that, for certain physical systems and for certain 
models and regimes of the coupling with the environment, one 
is able to estimate the typical time-scale within which the ef- 
fects of decoherence can be neglected. Hence one can deter- 
mine, in principle, a balanced working time for these systems. 
At this point one should actually check whether this balanced 
working time guarantees a suitable robustness of the quantum 
gate against parametric noise, namely, whether the effects of 
this kind of noise on the fidelity of the non-ideal quantum gate 
with respect to the ideal one can be neglected or not. 

Recently, a new ingredient has been proposed for the im- 
plementation of a holonomic quantum gate l24ll (see also l25l 
12610 where the authors have observed — for the model of a 
ion-trap quantum gate proposed by Duan et al. |27ll — the 
existence of a optimal working time, namely, of a specific op- 
erational time for which the non-ideal (i.e. finite-time) gate 
behaves exactly as the ideal (i.e. adiabatic) gate; they show, 
furthermore, that over the optimal working time the effects of 
the environment are negligible. Thus, such a optimal working 
time turns out to be also a balanced working time. 
We stress that, anyway, the fact that the non-ideal gate be- 
haves, in correspondence to the optimal working time, as the 
ideal one cannot be used to rule out the influence of parametric 
noise on the base of the standard geometric argument. Indeed, 
one should not expect that, perturbing the loop in the param- 
eter space, the non-ideal gate will still mimic the behavior of 
the ideal one. Hence one cannot apply, in principle, the stan- 
dard geometric argument to support the robustness of this kind 



of holonomic gate against parametric noise. A critical analy- 
sis of this simple, but somehow subtle, issue is the main aim 
of the present contribution. 

In conclusion, we think that the impact of parametric noise 
on holonomic gates is still an open problem and one is not 
legitimated, in general, to state the robustness of non-ideal 
holonomic gates against this kind of perturbations on the base 
a generic geometric argument. In our present contribution, 
we will try to illustrate this assertion by means of quantita- 
tive arguments, focusing on the ion-trap model proposed by 
Duan et al. 12711 . Even if other models have been proposed in 
the literature 12811 . the model of Duan et al. is probably the one 
most extensively studied also with reference to different phys- 
ical systems, as Josephson junctions lf29ll and semiconductor 
quantum dots 0"(ill . and can be regarded as a reference point 
for the subject. 

The paper is organized as follows. In Sec. [TT] the model 
Hamiltonian is introduced which will serve as a case study. In 
Sec. [HI] the behavior of the considered system in presence of 
several models of parametric noise is discussed. In Sec. HVIthe 
obtained results are analyzed and commented. Conclusions 
and remarks are presented in Sec. [V] 

II. A CASE STUDY 

As a case study, here we consider the single-qubit non 
Abelian gate that was proposed in j27il . The model under con- 
sideration can be physically realized as, for instance, a trapped 
ion with two degenerate ground (or metastable) states |0) and 
|1) which play the role of the computational basis. A quasi 
degenerate ancillary state \a) and an excited state |e) are also 
needed (the scheme is drawn in Fig. [Tf a ))- The low energy 
states are supposed to be independently coupled with the ex- 
cited state, such that the interaction picture Hamiltonian in the 
rotating frame reads as follows: 

H(r) = H(x, y, z) = |e) [x{0| + y(l| + z(a\] + h.c. . (1) 

The, in general complex, parameters x,y,z are related to three 
independent Rabi frequencies corresponding to, in general de- 
tuned, laser beams with different energies and polarization. In 
an ideal experiment, however, the laser beams are assumed to 
be resonant with the corresponding transitions and the param- 
eters are constrained to take values on a two-sphere. It is thus 
convenient to introduce polar coordinates: 

x = Q sin 1} cos tp 

y = rising sin tp . (2) 

z = cos d 

The spectrum of ([TJ) is threefold: a = {0, ±f£}, with the null 
eigenvalue which is doubly degenerate. The two degenerate 
eigenstates with vanishing energy can be chosen as follows: 

|-0o) = cos^ (coscp|0) + siny>|l)) — sint9|a), 
\tpi) = — sinyj|0) + cos <p\l). 

An analysis of the holonomy associated to the Hamiltonian 
<[TJ in correspondence with the doubly degenerate subspace 
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FIG. 1: Structure of the atomic levels and resonant lasers (a); unper- 
turbed loop © in the parameter manifold (b). 



shows that a closed path with starting point $ = corre- 
sponds to a non Abelian holonomy W = exp [— ia y uj), where 
<7 y = —i (|0)(1| — |1}(0|) is the Pauli matrix in the computa- 
tional space and co is the solid angle spanned by the param- 
eter 'd(s) and tp(s). This geometric character of the dynam- 
ics in the adiabatic limit is at the heart of the usual argument 
in favor of the robustness of holonomic quantum computa- 
tion. A stochastic noise in the control parameter can modify 
the details of the loop but, for a sufficiently great number of 
cycles of the noise during the system evolution, the fluctua- 
tions in the swept solid angle are consider to become negligi- 
ble (see lfl6l[r7ll and references therein). 

Here we consider the closed path in the parameter manifold 
that was studied in J3- For s G [0, 1] we take (see Fig. 02b)): 



<p(s) 



3stt/2 


s G [0,1/3] 


tt/2 


se [1/3,2/3] 


371-/2 (1 - s) 


s e [2/3,1] 





a G [0, 1/3] 


Stt/2 (,-§) 


s G [1/3,2/3] 


tt/2 


s G [2/3, 1] 



(4) 



The solid angle related to the loop (0J is ui = tt/2, hence 
the corresponding holonomic gate is W = —i<j y . As was 
observed in 12411 . the remarkable property of this path is that 
it presents perfect revivals of the gate fidelity at finite oper- 
ational time. The same behavior was predicted for all the 
loops constructed by moving from the north pole to the equa- 
tor through a meridian and back to the north pole through an- 
other meridian with piecewise constant velocity. In the case 
of the loop (01 there is a perfect revival of fidelity in corre- 
spondence of the operational times: 



k = S^iefc 2 -!, k = 1,2,. 
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(5) 



In the following we are mostly concerned with the first opti- 
mal operational time r* . 

To conclude this section we notice that a geometric phase 
appears in correspondence to a non adiabatic cyclic dynam- 
ics lfl4j[l5ll . In particular, for our case study, it happens that, 
in correspondence to an optimal operational time, the evolu- 
tion becomes cyclic and the acquired geometric phase is equal 
to the adiabatic holonomy. 



III. MODELS OF NOISE 

In order to study the robustness of non ideal holonomic 
gates, we consider the response of the system under paramet- 
ric noise in the ideal loop ©. In order to quantify the robust- 
ness of the gate, the noisy finite time evolution of the system 
is solved with analytic or numerical methods and the average 
gate fidelity is calculated. In the following, several models 
of noise are taken into account: in Sec. lIIlAl we consider the 
response of the system under a monochromatic perturbation 
of the three Rabi frequencies in (Q3; in Sec. MI Bl we consider 
a model of noise expressed by a random step function in the 
angular variables (O on the sphere; finally, in Sec. IHICI we 
discuss the response of the system under a random perturba- 
tion in the three Rabi frequencies. 



A. Monochromatic perturbation 

In this section we consider the behavior of the system in 
presence of a small perturbation in the parameters which can 
be viewed as a probe function used to test the stability of the 
gate, in particular we concentrate our attention to the case of 
a monochromatic probe function. A generic noisy path can be 
written as follows: 



r n (t) = r(t) + e(t), t€ [0,r], 



(6) 



where the vector r(t) describes the unperturbed loop and e(t) 
is a three component vector including the perturbation of the 
path. We have chosen a monochromatic perturbation at fre- 
quency T] and considered a noisy path obtained from d2j and 
©: 



C,t,0i) = x(s) + e „e'"' TS +* 
y n {s;r],e ni T,<f)2) = y(s) + > 



(7) 



^(s;r/,e^r,^ 3 ) = z(s) + e n e^ TS +^ 



where r n (s) = (x„(s), y n (s), z n (s)), 4> = (01,02,03) are 
random phases uniformly distributed in [0, 27r) and e v is the 
strength of the noise (chosen to be equal for the three compo- 
nent). Notice that this model of noise acts on both the ampli- 
tude and the de-tuning of the lasers. From (O it is clear that 
at finite operational time the perturbation does not reduces to 
a geometric perturbation of the loop in the parameters space 
since the perturbed path itself depends on the operational time. 
In presence of noise, different values of the operational time 
r correspond to different loops in the parameters manifold. 

For given values of rj, e v , t and </>, we consider the solution 
of the Schrodinger equation 

V^(s;r?,e^,0) = -iTH(r u (s))V T (s; <q, e n , 0), s G [0,1]. 

(8) 

where, in presence of noise, the re-scaled Hamiltonian 
H(r n (s)) depends on r too. Since we are mainly interested 
in the transformation emerging at the end of the loop, we set 
V T (v, e r,, 4>) = V T (l;r], e v , 0). 

Notice that, for all practical purposes, taking the average on 
the random phases corresponds to the action of the completely 



FIG. 2: Average gate fidelity as a function of the adimensional op- 
erational time Qt for several noise frequencies for the model in 
Sec. IIII Al Black boxes: e v — 0; circles: e v = 0.151, 77 = O.lfl; 
triangles: e v = 0.1Q,, r\ = 0.2S1; squares: e v — 0.1Q,, r\ — 0.3S1. 



FIG. 3: (Color on line.) Average gate fidelity at the first optimal 
operational time as a function of adimensional noise frequency (77/fi) 
and amplitude (e v /Q) for the model in Sec. IIII Al 



positive map 

P — S{p) 



(2tt) 3 



d(t>V T (f], e v ,(f))pV T (r),e n , </>) f . (9) 



This completely positive map has to be compared with the 
ideal adiabatic unitary dynamics, to do that, we have evaluated 
the average gate fidelity 



F 



dip(ip\W^£(\^)(ip\)W\ip) 



(10) 



where dijj indicates the normalized Fubini-Studi metric on 
pure states. This has been computed by means of the formula 
(see Ref. Q): 

F = i + ^tr [WW*£ (Po(0))]+-£ tr [WajWiSfa)] 

i=i 

(11) 

where <7j are the Pauli matrices in the computational subspace. 

For several values of rj, e n and <p, Eq. ^ is numerically 
solved using the relation: 



- iT H(r n (k/N))± 



V T ( V ,e v ,cl>)=lm L l[ k=oN exp 

(12) 

Where Y[ stands for the path ordered product. The effective 
completely positive map (0 is evaluated taking the average 
over 50 random choices of the phases <fi. Figure |2] shows the 
estimated gate fidelity ( fTTT i plotted as a function of the adi- 
mensional operational time fir, for several values of the noise 
amplitude and frequency. The unperturbed dynamics corre- 
sponds to en = and can be compared with the analytical 
results in 112411 . it exhibits perfect revivals of the average gate 
fidelity at finite time, in particular the first optimal operational 
time is Orj* ~ 18.25. The numerical results show that the pat- 
tern of gate fidelity as a function of the operational time can 
be completely different in presence of noise. 



The average gate fidelity at the first optimal operational 
time Ti in presence of parametric noise is plotted in Fig. [3] 
as a function of both amplitude and frequency of the noise. 
This plot suggests that the gate is indeed robust also for rather 
large noise amplitude (e v = 0.4S1). It is worth to notice that 
this is true unless the perturbation frequency is in a particular 
range approximatively about rj ~ 0.15S1. The presence of a 
typical frequency scale in the pattern of the fidelity is a fea- 
ture that will be reencountered in the other models of noise 
considered below. 

We have also studied, with the same methods, the response 
of the system in presence of analogous perturbations which 
have different symmetries. We have considered the case in 
which only the real part of (0 is taken; in this case the pertur- 
bation acts only in the amplitude of the coupling but not in the 
de-tuning. We have also analyzed the case of a perturbation 
which is square wave shaped; in this case a probe function is 
identified by its half period and initial phase. In both cases 
the corresponding patterns of the average gate fidelity are ex- 
actly analogous to the one shown in Fig. [3] This leads to the 
conclusion that the pattern of fidelity is largely independent of 
the details of the chosen probe function and a rather general 
behavior as function of the typical frequency is observed. 

Analogous results are also found for other loops of the same 
kind, such as the loop with the angle <p varying from to 7r/4 
in (O which is related to the Hadamard gate. 



B. Random noise on the sphere 

In this section we consider a model of noise which pre- 
serves the spectrum of the unperturbed Hamiltonian (fl3, be- 
cause of its symmetries an analytical solution of the noisy dy- 
namics is available. 

In ll24ll it was shown that the evolution operator can be eval- 
uated without approximation in several situations. Referring 
to the model in Eq. <£TJ, it is possible to evaluate the evolution 
operator along any segment on the parametric sphere as far as 
one of the parameters (z?, ip) is kept constant. In particular, re- 
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ferring to the case in Fig. 02b), the loop is composed by three 
segments and along each of them the previous condition is sat- 
isfied. Thus one can demonstrate Il24ll that the total evolution 
operator can be splitted in the form 

U(t) = U 3 (t 3 )U 2 (t 2 )U 1 (t 1 ), (13) 

where r is the total time evolution and 73, i = 1, 2, 3 are the 
times needed for covering each segment (for simplicity we 
suppose that the speed of the evolution are constant in each 
segment); moreover, the intermediate Ui's can be explicitly 
calculated ll24ll . Their form is very peculiar and it is possible 
to see that, in terms of the parameters ip), one can write 

f/i(*i) = E/i(tfi(*i),^i), (14) 
U 2 (t 2 ) = U 2 (-& 2 ,<p 2 (t 2 )), (15) 
U 3 (t 3 ) = U 3 (Mh)^ 3 ), (16) 

where < t{ < Ti and ip±, fl 2 and ip 3 are the constant values 
of the parameters during the evolution along the segment 1, 2 
and 3 respectively. 

We want to use these results for gaining information about 
the influence of the noise. We will therefore consider the fol- 
lowing model: every U{ is splitted in N evolution operators 
Uf evolving for a time r step (a sub-segment) such that 

N = T t /T step . (17) 

The evolution in the segment i reads 

N 

U i {T i ) = l[Ui(T s < ep ), (18) 

3=1 

In each sub-segment one of the sphere parameters is constant 
and the other evolves (we are moving on meridians or paral- 
lels). We add a random component to the constant parameter 
while the other is not affected. In other words, we are includ- 
ing a transverse component. We also suppose that the trans- 
verse evolution operator is equal to the identity (the "switch" 
is infinitely fast). This way we have splitted the evolution 
on a single meridian (parallel) in a sequence of evolutions 
of shorter meridians (parallels). Using Eq.s (fT4b-([T8l) we can 
write 

N 

U 1 (r 1 ) = l[U((Mt{)^i+e i ), d9) 

3=1 

N 

U 2 (t 2 ) = JJ Ui (02 + & , y 2 (4) ) , (20) 

3=1 

N 

Us(r 3 ) = l[ Ui(M4)><Ps (2D 

3 = 1 

where (j - l)r step < t \ < jT step and £f G [-7,7] are ran- 
dom variables uniformly distributed in the chosen interval 
(?' = 1, 2, 3). We stress again that each operator in the de- 
composition has a (large and not transparent) analytical ex- 
pression. Using this model we have computed the average 




FIG. 4: (Color on line.) Average gate fidelity at the first optimal 
operational time as a function of the re-scaled adimensional noise 
amplitude -y s = 7/(71-/2) and the adimensional parameter (f2T slep ) - 
for the model in Sec. lIIIBl 



gate fidelity at the first optimal operational time r* by means 
of Eq. ( fTTT i (and averaged over 50 realizations of the random 
process). The result is shown in Fig. |4j F is plotted as a func- 
tion of the noise amplitude 7 (re-scaled with the maximum 
value of the parameter for the loop in Fig. Q2b) i.e. 7r/2) and 
the parameter (ilr step ) _1 (characterizing the frequency of the 
noise); notice that, due to Eq. ( fTTI ), at fixed value of the op- 
erational time r, higher values of (Or step ) _1 correspond to 
a larger number of fluctuations. Also for this model the fi- 
delity exhibits a breakdown for small frequencies of the noise 
(which is in accordance with previous results). In particular, 
for (£lT step )~ < 0.5 the fidelity exhibits a minimum for any 
amplitude of the noise. Anyway, we notice that the deep of 
the fidelity is pronounced if the noise amplitude is one half 
the maximum value of the parameters; clearly, this situation 
corresponds to an unphysical scenario in which the control of 
the parameters is very poor. In all other situations the typi- 
cal values of F is very high. In the range of intermediate and 
large frequencies F quickly recovers the ideal behavior. 

It is interesting to compare the behavior of the gate at the 
first optimal operational time to the case of longer operational 
time in presence of noise, i.e., in the (approximated) adiabatic 
regime. It is possible to see I24I l25tl that the fidelity oscil- 
lations shown in Fig. [2] in absence of noise are strongly sup- 
pressed if k > 3 in Eq. (0 (we are near the adiabatic regime). 
A good approximation of the adiabatic regime can be already 
obtained for the fourth optimal operational time. Therefore, 
we have computed the average gate fidelity for firf ~ 75.21 
[(£lT step )~ ranges as in Fig. 2). The result is shown in Fig. [5] 
and can be directly compared to the plot in Fig. [4] First of all 
it is important to stress that the total number of fluctuations 
A^4 for J7t| ~ 75.21 is larger when compared to the number 
Ni for the first optimal working point f2rf ~ 18.25. From 
Eq. (O and Eq. (T% we have ~ 4.12^1. In apparent 
contrast to the intuition related to the usual argument of ro- 
bustness of holonomic gates we notice that, in the same range 
of frequencies of the non adiabatic case (and, therefore, for a 
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FIG. 5: (Color on line.) Average gate fidelity at the fourth optimal 
operational time as a function of the re-scaled noise amplitude 7 S = 
7/(tt/2) and the adimensional parameter (r2r st ep) _1 for the model 
in Sec. HlTBl 
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FIG. 6: Average gate fidelity as a function of the noise typical fre- 
quency for the noise model in Sec. IIII CI for the first four optimal 
operational times. Triangles, circles, full triangles and squares corre- 
spond respectively to the first, second, third and fourth optimal oper- 
ational time, e = 0.1S1. 



larger number of fluctuations), F reaches lower values. More- 
over, the adiabatic gate needs higher values of the frequency 
of noise for recovering the ideal behavior. We conclude that 
the (approximately) adiabatic (purely geometric) NOT trans- 
formation is more sensitive to parametric noise than the non 
adiabatic one. 



C. Random noise 

In this section we consider a model for a random perturba- 
tion of the loop which is not constrained to preserve the sphere 
in the parameter space. Taking in consideration the ideal loop 
here we study the noisy paths of the following kind: 

Xn{s;T step ,e,T) = X(s) +£l(s,T step ,T) 

yn(s;T step ,e,r) = y(s) + £ 2 (s, r step , t) , (22) 

Zn(s;T ste p,e,r) = z(s) +£3(s,T step ,T) 

where £i(s, T ste p, r) G [— e, e] are three real random variables, 
uniformly distributed in the chosen interval, which are piece- 
wise constant for (j — l)r step < sr < jr step . 

In order to study the behavior of the gate at finite opera- 
tional time, we have evaluated the average gate fidelity for a 
fixed value of the noise amplitude e = 0.1f2 as a function of 
the noise typical frequency (r2r s t C p) _1 in correspondence of 
the first four optimal working times. The results are shown in 
Fig-El T ne data plotted in this figure lead us to two considera- 
tions: first of all we notice again the unexpected result that the 
non-adiabatic optimal working times (the first, for instance) 
appears to be more robust than longer operational times (the 
forth optimal operational time, for instance); secondly, we ob- 
serve the same qualitative behavior of the pattern of fidelity 
for all the optimal operational times under study, this suggests 
the presence of a common mechanism which account for the 
cancellation of the effects of the noise. 

We have also analyzed the case of a noise which in- 
clude de-tuning by considering complex random variables 



£i( s ! T s tep, t). The result are completely analogous and the in- 
troduction of a noise in the de-tuning does not introduce new 
elements in the pattern of fidelity. 



IV. ANALYSIS OF THE RESULTS 

The aim of this section is to look for a physical explana- 
tion of the observed robustness of the considered finite time 
non-adiabatic gate. Due to the fact that all the models of noise 
induce the same qualitative behavior of the fidelity, in the fol- 
lowing we are going to consider in more details the model pre- 
sented in Sec. IIII Cl It is worth to notice that only for the first 
model the noise affects both the amplitude and the phase of the 
control parameters, while the other two models of noise con- 
cern only their amplitude. Nevertheless, it is a result of lfl7ll 
that the main contribution in the noisy dynamics is due to the 
component in the amplitude. 

As already recalled, the relevant parameter for the geomet- 
rical cancellation usually related to holonomic gates in the adi- 
abatic regime is the number of fluctuations of the noise during 
the gate operational time (denoted N). This effect is related 
only to the swept solid angle and is independent of the chosen 
operational time. If the number of cycles of the noise is large 
enough, the fluctuations in the solid angle spanned by the loop 
are expected to become negligible. To be more specific sup- 
pose that, after a noisy loop, the swept solid angle is u> and the 
mean square over the realizations of the noise is (Auj 2 ). In 
Fig- ID the mean square is plotted as a function of the number 
of cycles of the noise; since, in the adiabatic limit, the gate 
depends only on the swept solid angle, the fluctuations of the 
gate are expected to have the same behavior as the fluctuations 
in the solid angle. 

As already explained in the previous section, Fig. [6] shows 
the average gate fidelity as a function of the adimensional typ- 
ical noise frequency (f2r s t cp ) _1 for several values of the evo- 
lution time which correspond to the first four optimal oper- 
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FIG. 7: Fluctuations in the solid angle spanned by a noisy loop as a 
function of the number of perturbations of the noise N, for the noise 
model in Sec. lmCl 6 = O.lfi. 
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FIG. 8: Average gate fidelity as a function of the number of fluctu- 
ations of the noise TV for the noise model in Sec. ImCl for the first 
four optimal operational times. Triangles, circles, full triangles and 
squares correspond respectively to the first, second, third and fourth 
optimal operational time, e = O.lfi. Compare with Fig.[6]and|7] 



ational times. The plot shows an analogous behavior of the 
fidelity as a function of the typical noise frequency indepen- 
dently of the particular value of the operational time; more- 
over, the minimum of the fidelity is reached in correspondence 
of (OT s tep) -1 — 0.5 for all the values of the operational time 
considered. In order to cast some light on the nature of the 
cancellation effect, the same data are plotted in Fig.[8]as func- 
tions of the number of fluctuations of the noise (notice that 
N = f2r(ilr s tc P ) _1 ). A direct comparison of figures [6] and 
|8] suggests that the relevant quantity which accounts for the 
mechanism of cancellation of the effects of the noise is its 
typical frequency (rir s t C p) _1 and not only the number of fluc- 
tuations N. On the other hand, the fluctuations of the solid 
angle around the ideal value (tt/2) start to be negligible for 
N > 20; a comparison with the curve for the fourth optimal 
working point (squares in Fig. [8} suggests that the recovery 
of the fidelity for long cyclic evolution times is given also by 
geometric cancellation. 

For non adiabatic evolution times one can imagine the ex- 



istence of a different mechanism which accounts for the ob- 
served cancellation of the noise effects for sufficiently fast 
noise which is related to a dynamical instead of geometrical 
cancellation. A dynamical effect could not be directly related 
to the swept solid angle: in this case the relevant parameter is 
expected to be the typical time of the noise r s te P and a dynam- 
ical cancellation of the noise should appear if its typical fre- 
quency is sufficiently large compared to the system frequency, 
namely (rir s t ep ) _1 3> 1. Of course this condition implies, 
for fixed operational time r, that N 3> 1 (the usual condition 
for geometric cancellation); nevertheless, as Fig. [6] shows, a 
cancellation of the noise effects appears on a frequency scale 
(fiTstop) -1 ~ 1 independently of the chosen value of the op- 
erational time, thus suggesting a dynamical mechanism for the 
noise cancellation at least for the first four optimal operational 
times. 

The fact that in the non adiabatic regime the robustness has 
a dynamical origin can also explain why the minimum value 
of the fidelity tends to decrease for increasing values of r*: 
if the geometric cancellation is not present, the noise is less 
effective in disturbing the system when the evolution time is 
short. 



V. CONCLUSIONS 

In this paper we have considered the influence of paramet- 
ric noise on the efficiency of a non adiabatic holonomic gate 
which is known to be robust in the ideal case. Three models of 
parametric noise have been discussed in the case of finite op- 
erational time. The average gate fidelities for all the models of 
noise considered here present an analogous qualitative behav- 
ior. For each of the three models the non ideal gate presents 
a breakdown of the average gate fidelity for small frequencies 
of the noise (compared to the system Bohr frequency), while a 
high value of the fidelity is reached for noise with higher fre- 
quencies. This can lead to say that the presence of a "resonant 
frequency" for the breakdown of F is a general feature of any 
model of parametric noise. 

We want to stress again that the usual argument in favor 
of the robustness of holonomic quantum computation is based 
on the purely geometric nature of the holonomy group that de- 
scribes the adiabatic transformations. Since the dynamics has 
a completely geometric character only in the adiabatic limit, 
the robustness of adiabatic gates is, in this sense, just a con- 
sequence of the adiabatic theorem. Despite these considera- 
tions, our calculations show that, at least in certain situations, 
the first optimal operational time can be preferable to longer 
operational times with regards to the robustness of the corre- 
sponding gate against parametric noise. 

Nevertheless, our results lead to the conclusion that the ob- 
served revivals of the fidelity for sufficiently fast noises is 
mainly due to dynamical instead of geometrical effects. Our 
conclusion is that, in the range of operational times consid- 
ered here, the observed cancellation effects are mainly re- 
lated to a dynamical average over fast oscillations of the noise 
(f^Tstcp) -1 3> 1 and there is no relevant connection with the 
swept solid angle which plays a crucial role for the usual ar- 
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gument in favor of robustness of the holonomic computation 
in the adiabatic regime. 
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